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ABSTRACT 

Natural,  mixed  and  forced  convection  liquid  cooling  of 
discrete  heat  sources  in  a  vertical  channel  were  investigated 
experimentally.  Ten  heat  sources  were  flush  mounted  to  one 
wall  of  a  water  filled  channel,  while  the  opposite  wall  was 
unheated.  Measurements  of  heater  surface  temperatures  were 
made  for  a  channel  Reynolds  number  range  of  0  -  17  00  and  heat 
flux  range  of  1050  -  4500  W/m2.  In  the  mixed  convection 
regime,  both  aiding  and  opposed  flow  were  examined. 
Temperature  patterns  on  the  heated  surface  were  visualized 
using  liquid  crystals. 
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I .  INTRODUCTION 

A.   ELECTRONIC  COOLING 

With  the  invention  of  the  transistor  in  1947,  the  present 
computer  age  was  born.  The  limit  on  computer  design  was 
circuit  technology.  As  this  technology  grew  and  engineers 
developed  ways  to  pack  more  devices  on  each  chip,  heat 
dissipation  became  a  limiting  design  factor. 

Efficient  heat  dissipation  from  highly  packed  circuits  is 
vital  to  reliable  chip  operation.  For  every  10°C  rise  in 
junction  temperature  above  85°C,  reliability  suffers  by  a 
factor  of  two  [Ref  1:  p.  35].  The  importance  of  adeguate 
cooling  is  better  visualized  by  considering  a  typical  5  mm  by 
5  mm  chip  dissipating  10  W  of  power.  This  gives  an  average 
heat  flux  of  4xl05  W/m2 — just  two  orders  of  magnitude  less  than 
that  on  the  sun's  surface. 

Electronic  cooling  is  evolving  with  advances  in  semi- 
conductor fabrication  technology.  As  Small  Scale  Integration 
(SSI)  has  developed  into  Ultra  Large  Scale  Integration  (ULSI) , 
so  has  air  cooling  technology  developed  into  conduction  and 
direct  liquid  cooling. 

1.  Air  Cooling  Technology 

For  low  power  densities,  air  cooling  is  commonly  used 
to  cool  electronic  components.  Within  air  cooling  technology, 


the  simplicity  of  natural  convection  has  given  way  to 
sophisticated  forced  convection  designs  in  order  to  dissipate 
the  heat  from  high  power  systems. 

In  order  to  enhance  heat  removal  capabilities, 
designers  have  had  to  rethink  module  and  system  designs.  The 
Hitachi  SiC  RAM  and  IBM  4301  modules  are  two  examples  [Ref  2: 
p.  294].  These  designs  have  attempted  to  minimize  internal 
and  external  heat  flow  resistances  across  the  module  as  well 
as  made  use  of  parallel  flow  schemes  to  reduce  the  temperature 
rise  across  the  system.  In  addition  to  the  above,  other 
enhancement  techniques  being  explored  and  used  include: 
impinging  air  flow  schemes,  turbulators,  and  water-cooled 
"interboard"  heat  exchangers  [Ref  2:  pp.  293-296]. 

While  much  has  been  done  to  enhance  air  cooling 
technology,  these  developments  have  failed  to  keep  pace  with 
the  heat  dissipation  requirements  of  modern  electronic 
equipment. 

2.  Conduction  Cooling  Technology 

IBM's  Thermal  Conduction  Module  (TCM)  is  an  example  of 
keeping  pace  with  today's  higher  power  densities.  The  key 
elements  in  this  technology  include:  a  modular  cold  plate 
subsystem,  a  coolant  distribution  loop,  heat  rejection  loop 
and  appropriate  thermal  sensing  and  control  hardware  [Ref  2: 
p.  298].  Briefly,  the  heat  path  goes  from  the  chip  surface 
through  a  spring  loaded  piston  and  housing  to  the  cold  plate. 


Enhancements  to  this  technology  center  around  reducing  the 
thermal  resistance  and  include:  optimum  geometry  of  the 
piston,  housing  and  cold  plate;  improved  thermal  fillers 
across  physical  gaps;  and  optimum  flow  rates.  [Ref  2:  pp.  298- 
299]  . 

A  promising  enhancement  technique  which  can  be 
classified  as  either  conduction  or  direct  cooling  is  the 
microgroove/microchannel  technique.  First  proposed  by 
Tuckerman  and  Pease  [Ref  3:  pp.  126-129],  this  technique  would 
significantly  reduce  the  conductive  resistances.  The  cold 
plate,  in  effect,  is  brought  to  the  under  side  of  the 
integrated  circuit  chip. 

3 .  Direct  Liquid  Cooling 

Because  of  higher  density  and  specific  heat  of  liquids 
compared  to  gases,  liquid  immersion  cooling  has  a  clear  heat 
transfer  rate  advantage  over  air  cooling  technology. 
Immersion  and  circulation  of  coolant  through  microgrooves  in 
the  chip  eliminate  the  interfacial  resistance  that  restricts 
conduction  cooling  technology. 

The  Cray-2  and  ETA-10  supercomputers  are  two  examples 
of  the  use  of  direct  liquid  cooling.  With  the  success  of 
these  systems,  researchers  and  designers  have  had  an  increased 
interest  in  this  method  of  cooling.  New  fluids  are  being 
developed  to  meet  the  strict  chemical,  electrical  and  thermal 
requirements.   Additionally,  as  with  any  new  technology,  a 


knowledge  base  of  data  is  beginning  to  be  developed.  Although 
much  research  has  been  conducted  for  air  cooling,  relatively 
less  has  been  done  with  direct  liquid  cooling.  Due  to  the 
large  differences  in  the  thermophysical  properties,  the  air 
cooling  data  in  the  literature  is  not  directly  applicable  to 
liquid  cooling.  The  modes  of  direct  liquid  cooling  are: 
natural,  forced  and  mixed  convection;  phase  change;  falling 
film;  and  jet  impingement. 

B.   PREVIOUS  RESEARCH  ON  SINGLE  PHASE  LIQUID  COOLING 

In  this  section  and  the  one  following,  a  brief  review  of 
single  phase  liquid  cooling  for  natural,  forced  and  mixed 
convection  is  presented.  In  order  to  establish  a  knowledge 
base,  it  is  conventional  to  start  with  simple  configurations, 
attempt  to  fully  understand  the  behavior  observed  and  then 
expand  this  understanding  to  more  complicated  situations  by 
both  computer  simulation  and  further  experimentation.  For 
natural,  mixed  and  forced  convection,  typical  configurations 
simulating  electronic  cooling  applications  include  channels, 
isolated  surface  and  full  or  partial  enclosures. 

Park  and  Bergles  have  investigated  natural  convection  for 
discrete  heat  sources  on  isolated  surfaces  [Ref  4:  pp.  90-96]. 
By  varying  heater  size  and  arrangements  they  found  that 
shorter,  wider  heaters  provided  higher  heat  transfer 
coefficients;  protruding  heaters  performed  better  than  flush 
mounted  heaters;  downstream  heater  heat  transfer  coefficients 


were  lower  than  upstream  ones;  and  staggered  vice  in  line 
configuration  was  better  for  heat  transfer. 

Joshi  and  Kelleher  [Ref  1:  pp.  35-39]  have  described  a 
number  of  natural  convection  liquid  immersion  cooling  studies 
for  the  three  configurations  identified  above.  Both 
experimental  and  computational  studies  are  described. 
Applications  of  natural  convection  cooling  of  actual  packages 
are  addressed. 

Keyhani  et  al.  [Ref  5:  pp.  616-623]  considered  two 
dimensional  natural  convection  from  eleven  discrete  flush 
mounted  heat  sources  on  a  vertical  wall  in  a  rectangular 
enclosure.  Heat  transfer  data  and  flow  visualization  revealed 
primary,  secondary  and  tertiary  flow  patterns  and  indicated 
that  stratification  was  an  important  factor  influencing  the 
temperature  of  the  heated  sections.  These  results  also 
revealed  higher  Nusselt  numbers  for  discrete  heating  than  for 
a  fully  heated  vertical  wall. 

Forced  convection  liquid  cooling  research  has  centered 
around  the  earlier  work  done  by  Tuckerman  and  Pease  on 
microgroove/microchannel  configurations.  Kishimoto  and  Sasaki 
[Ref  6:  p.  358]  proposed  methods  to  reduce  temperature 
variations  across  silicon  chips  by  using  diamond  shaped 
channels.  Phillips  et  al.  [Ref  7:  pp.  227-248]  provide  models 
for  heat  transfer  and  pressure  drops  for  various  aspect  ratio 
channels  for  laminar  and  turbulent  flow  regimes.   Phillips 


[Ref  8:  pp.  31-48]  also  presented  a  thermal  resistance  model 
for  predicting  microchannel  cooling  performance. 

C.   MIXED  CONVECTION 

In  many  applications,  both  natural  and  forced  convection 
are  important.  Mixed  convection  deals  with  the  flow  regime 
where  characteristic  buoyant  and  forced  convection  velocities 
are  comparable.  The  effects  of  buoyancy  forces  and  the 
specific  parametric  ranges  where  they  become  important  in  the 
overall  heat  transfer  process  are  highly  dependent  on  the 
geometry  and  boundary  conditions. 

Mixed  convection  combines  the  natural  and  forced  flow 
effects  in  the  form  Gr/Ren,  where  the  Grashof  number  (Gr)  is 
a  measure  of  the  ratio  of  buoyancy  to  viscous  forces  and 
Reynolds  number  (Re)  indicates  the  ratio  of  inertia  to  viscous 
forces.  The  exponent  'n'  is  a  function  of  geometry,  fluid, 
and  heating  condition.  As  this  nondimensional  ratio 
approaches  zero,  forced  convection  dominates;  as  it  approaches 
infinity,  natural  convection  dominates.  Mixed  convection  may 
be  aiding  or  opposed.  Aiding  implies  free  stream  flow  in  the 
direction  of  buoyancy  forces;  opposed  convection  implies  free 
stream  flow  opposite  to  buoyancy  forces. 

Cheng  et  al.  [Ref  9:  p.  375]  obtained  analytical  solutions 
for  buoyancy-assisted,  fully  developed  laminar  mixed 
convection  flow  in  a  channel  for  three  combinations  of  thermal 
boundary  conditions — isof lux-isof lux,  isof lux-isothermal  and 


isothermal-isothermal.  For  these  conditions,  effects  of  flow 
reversal  and  heat  transfer  characteristics  were  studied. 

Maheny,  Incropera  and  Ramadhyani  studied,  experimentally 
and  numerically,  mixed  convection  from  an  in-line,  four  row 
array  of  twelve  discrete,  square,  flush  mounted  heat  sources 
in  a  liquid  filled  horizontal  channel.  They  determined  the 
variations  of  row-average  Nusselt  vs  Reynolds  number  and  found 
that  heat  transfer  may  be  enhanced  due  to  buoyancy  effects 
thereby  reducing  the  required  pumping  power.  [Ref  10:  pp. 
1233-1243] 

Mori  and  Ohbuchi  performed  experiments  on  the  dependence 
of  heat  transfer  on  heat  flux  for  downward  flows  in  a  vertical 
straight  circular  tube  for  Reynolds  numbers  less  than  2  000. 
Their  study  showed  that  for  small  heat  fluxes  the  temperature 
and  velocity  gradients  at  the  tube  wall  were  decreased  by 
buoyancy  forces  and  that  heat  transfer  was  reduced.  At  higher 
heat  fluxes,  buoyancy  and  velocity  fluctuations  appeared  at  Re 
>  500  resulting  in  enhancemented  heat  transfer.  [Ref  11:  pp. 
1231-1237] 

D.   OBJECTIVES  OF  PRESENT  STUDY 

This  study  experimentally  investigated  natural,  mixed  and 
forced  convection  heat  transfer  in  a  water  filled  vertical 
rectangular  channel  with  discrete  heater  strips  on  one  wall. 
Specific  objectives  were: 


•  To  investigate  natural  convection  effects  and  compare  the 
heat  transfer  data  with  known  correlations. 

•  To  investigate  heat  transfer  characteristics  for 
opposed  mixed  convection. 

•  To  investigate  heat  transfer  characteristics  for 
aiding  mixed  convection. 


II.   EXPERIMENTAL  APPARATUS 

A.   TEST  CHANNEL  ASSEMBLY 

The  test  channel  assembly  is  similar  to  the  one  used  by 
Syring.  [Ref  12]  It  consists  of  three  main  systems:  fluid 
circulation;  power  distribution;  and  data  acquisition  (Figure 

1). 

1.   Fluid  Circulation  System 

Water  is  circulated  in  a  closed  system  consisting  of: 
a  vertical  rectangular  channel  and  associated  tubing;  positive 
displacement  pump;  and  flow  meter.  The  vertical  channel 
consists  of  two  1  cm  thick  plexiglass  boards  2.0  m  by  0.25  m 
separated  by  a  1  cm  thick  plexiglass  spacer  with  a  2.5  cm  wide 
and  1  mm  thick  rubber  gasket.  At  the  top  and  bottom  of  the 
channel  are  surge  tanks.  Within  the  channel  are  the  test 
section  and  two  flow  straightening  regions.  The  top  flow 
straightening  region  (encountered  by  the  flow  prior  to 
reaching  the  test  section  during  opposed  mixed  convection)  is 
located  1.0  m  above  the  top  of  the  test  section.  The  other 
region  used  for  conditioning  the  flow  during  the  aiding 
convection  studies  is  located  0.3  m  below  the  bottom  of  the 
test  section. 

The  test  section  constitutes  a  20  cm  x  19  cm  portion 
of  one  wall  of  the  vertical  channel.   Ten  0.5  cm  x  20  cm  and 
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Figure    1.    Overall    schematic   of    experimental    apparatus 


10 


0.120  cm  deep  grooves  were  milled  horizontally  across  the  test 
surface  at  spacings  of  1  cm.  These  grooves  served  as 
thermocouple  beds  for  each  heating  strip.  A  second  set  of 
ten,  wider  but  shallower,  grooves  (1.0  cm  x  2  0  cm  and  0.08  cm 
deep)  were  milled  on  top  of  each  of  the  thermocouple  grooves. 
These  grooves  housed  the  flush  mounted  heating  strips.  An 
example  of  a  typical  heating  groove  surface  is  shown  in  Figure 

2.  Finally,  a  vertical  1.4  cm  x  29  cm  and  0.120  cm  deep 
groove  was  milled  to  house  the  electrical  leads  from  the  ten 
heater  strips.  A  front  view  of  the  test  surface  is  shown  in 
Figure  3.  A  side  view  of  the  channel  is  shown  in  Figure  4. 
Before  mounting  the  heater  strips  to  the  plexiglass,  five 
copper-constantan  (Cu-Co)  thermocouples  were  mounted  in  the 
deepest  groove  for  each  heater.  This  resulted  in  a  total  of 
fifty  thermocouples  for  the  entire  test  section. 
Thermocouples  for  each  heater  were  spaced  as  shown  in  Figure 

3 .  The  two  end  thermocouples  for  each  heater  were  made  with 
0.2  54  mm  diameter  wire;  the  middle  three  were  made  from  0.12  7 
mm  diameter  wire. 

All  of  the  thermocouples  were  mounted  using  a  low 
thermal  conductivity  fast  curing  exopy  (Omega  bond  100)  .  This 
served  to  effectively  reduce  the  conduction  losses  across  the 
plexiglass.  A  very  thin  layer  of  a  high  thermal  conductivity 
epoxy  (Omega  bond  101)  was  placed  in  the  area  of  the  heating 
strips  where  the  thermocouples  would  be  measuring.  The 
heating  strips  were  then  carefully  mounted  in  their  respective 
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Figure  2.  Typical  heating  surface  groove.   Dimensions  in 
cm. 
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Figure  3.  Front  view  of  test  surface.   Dimensions  in  cm 
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grooves  using  a  very  thin  coat  of  Omega  bond  100  to  ensure  a 
flush  mounted  heating  surface. 

A  Cole-Parmer  positive  displacement  gear  type  pump 
driven  by  a  variable  speed  motor  produced  Reynolds  numbers  in 
the  range  of  interest  (0  to  1700)  .  An  Omega  turbine  flow 
meter  was  used  to  monitor  the  fluid  flow  rates.  The  flow 
meter  was  calibrated  by  comparing  its  output  to  measurements 
of  guantities  of  water  collected  over  specified  times.  The 
maximum  uncertainty  in  flow  rate  was  0.5%.  Teflon  fittings 
and  tubing  provided  the  connections  between  the  vertical 
channel,  pump  and  flow  meter. 

2 .   Power  Distribution  System 

A  schematic  of  the  power  distribution  panel  is  shown 
in  Figure  5.  Power  to  the  individual  heaters  was  supplied  by 
a  Hewlett  Packard  HP  6269B,  0-40  VDC,  0-50  A  power  supply. 
This  supply  was  connected  to  a  common  bus  which  provided  power 
to  each  of  the  ten  heating  strips.  Each  heater  was  connected 
in  series  with  a  precision  resistor  (R  =  2.01  ±  0.005  Q) .  The 
data  acquisition  system  obtained  voltage  readings  across  the 
entire  bus  bar  (VT  )  and  across  each  heater  (VHeat)  .  From  these 
two  measured  values,  the  input  power  was  computed  using  the 
following  relationship: 
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3 .   Data  Acquisition  Assembly 

The  data  acquisition  system  consisted  of  a  HP  9000 
computer  and  a  HP  3852  data  acquisition  unit.  The  computer 
instructed  the  data  acquisition  system  to  monitor  temperatures 
from  desired  thermocouples  and  measure  the  power  delivered  to 
each  heater  strip.  Only  selected  thermocouples  were  scanned 
for  each  data  run,  depending  upon  the  information  desired. 
The  thermocouples  were  scanned  one  at  a  time  at  a  sampling 
rate  of  5  samples/second. 

B.   LIQUID  CRYSTALS 

The  surface  temperature  patterns  were  visualized  using  a 
0.762  mm  thick,  20  cm  X  20  cm,  Hallcrest  thermochromic  liquid 
crystal  (TLC)  sheet  bonded  to  the  test  surface  as  seen  in 
Figure  4.  The  crystals  on  this  sheet  displayed  shades  of 
green  and  blue  in  reaction  to  temperature  changes  in  the  range 
of  3  0  to  3  5°C.  The  observed  color  bands  were  compared  with 
the  average  inlet  temperature  TIN.  The  first  appearance  of 
green  and  blue  corresponded  to  3  0.5  °C  and  3  6.0  °C 
respectively. 
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C.   DISCUSSION  OF  SPECIAL  FEATURES 

Before  explaining  the  experimental  procedure  in  detail,  a 
general  discussion  of  certain  design  and  data  reduction 
procedures  is  presented. 

With  a  fixed  vertical  rectangular  geometry,  flow  rates  and 
power  inputs  were  varied.  Mean  flow  velocities  ranged  from  0 
to  0.08  m/s  and  heat  fluxes  from  1050  W/m2  to  4500  W/m2.  With 
these  varying  conditions,  experimental  data  runs  were 
performed  for  no  flow,  opposed  and  aiding  flows. 

For  laminar  conditions,  the  axial  length  for  fully 
developed  flow,  1E,  is: 


1E  =   0 . 16*a* 


'p.-x 


where,  Um  is  the  mean  velocity  and  'a'  the  channel  half  width 
[Ref  13:  pp.  185-186].  This  criterion  is  for  flow  between 
flat  parallel  plates.  Table  I  provides  sample  entry  length 
values  for  various  Reynolds  numbers.  Based  on  these  lengths, 
flow  straighteners  were  positioned  in  the  test  channel  for 
both  the  aiding  and  opposed  mixed  convection  studies. 

Free  stream  temperature  (TIN)  was  evaluated  by  taking  the 
average  of  two  Cu-Co  thermocouples  installed  on  the  inlet  side 
of  each  set  of  flow  straighteners.   Water  properties  were 
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TABLE  I.  CHANNEL  ENTRY  LENGTH  VALUES  FOR   VARIOUS 
REYNOLDS  NUMBERS 


Re 

GPM 

Um  (m/s) 

1E  (m) 

500 

0.7516 

0.0237 

0. 1052 

1000 

1.5640 

0.0493 

0.2105 

1500 

2.2963 

0.0724 

0. 3155 

1700 

2.6037 

0.082 

0.378 

evaluated  based  on  TIN.  The  center  temperature  (THTR)  was  used 
for  all  heater  measurements  because  heater  strip  temperature 
was  symmetrical  about  the  center  and  showed  only  small 
spanwise  variation  (within  a  0.2°C  band). 

Conduction  loss  through  the  substrate  to  the  outside 
environment  was  neglected.  For  a  heat  source  where  the 
surface  area-to-perimeter  is  large,  Baker  [Ref  14:  p.  172] 
suggests  that  the  substrate  conductivity  would  have  a  small 
effect  on  total  heat  transfer. 

D.   EXPERIMENTAL  PROCEDURE 

For  natural  convection,  data  were  taken  with  only  heater 
strip  1  energized  for  nine  different  heat  fluxes  (500,  1000, 
1500,  2000,  2500,  3000,  3500,  4000  and  4500  W/m2)  .  After  a 
steady  state  condition  was  reached,  five  minutes  of  data, 
sampled  at  a  sampling  freguency  of  5  hz,  were  taken  on  the 
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center  thermocouple  for  each  flux  level.  This  procedure  was 
repeated  for  heater  strips  2  and  3.  A  total  of  27  data  runs 
were  performed. 

For  opposed  and  aiding  mixed  convection,  all  heater  strips 
were  energized  at  four  different  power  levels  (1050,  2500, 
3500  and  4500  W/m2)  and  11  different  Reynolds  numbers  were 
considered  (100,  200,  300,  400,  500,  600,  800,  1000,  1200, 
1500  and  1700) .  After  a  nominal  steady  state  condition  for 
flux  and  flow  rate  was  achieved  (approximately  five  minutes) , 
a  five  minute  continuous  data  run  was  performed  on  all  center 
thermocouples.  For  a  given  flux  level,  data  runs  were 
conducted  in  a  high  to  low  Reynolds  numbers  seguence  (1700 
through  100)  .  A  total  of  44  data  runs  were  conducted  for 
opposed  and  44  data  runs  for  aiding  flow. 
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III.  RESULTS 

A.   DATA  ANALYSIS 

Neglecting  substrate  conduction  losses,  a  simple  energy 
balance  for  each  heater  strip  yields: 

QcONV    '       ^Heat 

where, 


®Heat 


V        -   V 

v  Tot  v Heat 


R 


*  v 

v Heat 


and 

Qconv  =  hMT 
In  the  above  equations,  h  is  the  convective  heat  transfer 

coefficient,   As  is  the  heater  surface  area  and  AT  the 

difference  between  THTR  and  TIN. 

The  average  Nusselt  number  (Nu)  was  calculated  as: 

at    h  *  L 
Nu  =   


Reynolds  number  (Re)  was  calculated  by 


Um    *    DH 

Re  =  -JH « 
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where,  DH  is  the  hydraulic  diameter  defined  as 

4  *  Ar 

H  p 


and  Ac  is  the  channel  cross  sectional  area  and  P  the  channel 
perimeter. 

Flux  based  Grashof  number  (Gr*)  was  calculated  by: 

KV2 


where,  P  is  the  volumetric  expansion  coefficient  for  the  fluid 
and  q"  the  heat  flux. 

Nondimensional  temperature  was  calculated  using: 

q  _  At*k 

q"*L 

The  uncertainties  for  Nu,  Re  and  Gr  are  5.2,  6.5  and  5.1% 
respectfully.  Sample  calculations  are  found  in  Appendix  A  and 
B. 

Data  reduction  for  this  study  resulted  in  the  following 
plots: 

•  Natural  convection:        Nu  vs  Gr*. 

•  Opposed  mixed  convection:   AT   vs   Re; 

Nu  vs   Re;  and 
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@   vs  Re. 
•   Aiding  mixed  convection:    AT   vs   Re; 

Nu  vs   Re . 

B.   NATURAL  CONVECTION  RESULTS 

The  primary  reason  for  looking  at  natural  convection  was 
to  establish  a  level  of  confidence  in  the  test  apparatus  and 
to  compare  present  data  with  two  published  constant-heat-flux 
vertical  wall  correlations. 

Figure  6  shows  a  comparison  of  the  present  data  with  Fujii 
and  Fujii  [Ref  15:  p.  320]  correlation  for  a  fully  heated 
vertical  plate. 

Nu    =   2^ —  =  f IL ]  '  '  *  (Gr\Pr)0-2 


This  correlation  for  the  local  Nusselt  number  applies  for 
laminar  flow  along  semi-infinite  vertical  flat  plate 
maintained  at  a  uniform  heat  flux.  Also  shown  in  Figure  6  is 
an  experimental  correlation  from  Joshi  and  Knight  [Ref  16], 
for  a  single  column  of  in  line  flush  mounted  heat  sources  on 
a  vertical  substrate  given  by  Nu=l.  11*  (Gr*)  °*12.  The  comparison 
is  shown  for  heater  strips  1,  2  and  3  of  the  present  study. 
The  local  Nusselt  numbers  were  evaluated  at  the  heater  strip 
vertical  midpoints  (0.005  mm). 
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The  heat  transfer  coefficients  for  the  heater  strips  were 
57%  below  the  Fujii  and  Fujii  and  36%  below  the  Joshi  and 
Knight  correlations.  The  power  law  exponents  (n)  for  Nux  = 
Constant*Grx*n  for  References  [15],  [16]  and  this  study  are 
0.2,  0.12  and  0.057  respectively. 

The  difference  between  Nux  for  References  [15]  and  [16] 
are  attributed  to  heater  size,  geometry  and  boundary 
conditions.  In  particular,  the  use  of  discrete  heat  sources 
in  the  latter  and  the  fact  that  a  uniform  flux  condition  is 
not  easily  achieved  with  small  heat  sources  are  likely  to  be 
the  primary  reasons  for  these  deviations.  The  present 
configuration  differs  from  Reference  [16]  in  that  the  heaters 
are  much  wider  and  also  the  opposite  channel  wall  is  likely  to 
have  an  influence.  The  self -consistency  of  the  present 
results  is  indicated  by  the  good  agreement  at  the  three 
locations — HS1,  2  and  3. 

C.   OPPOSED  MIXED  CONVECTION  RESULTS 

Four  different  heat  flux  settings  (1050,  2500,  3500  and 
4500  W/m2)  and  eleven  Reynolds  numbers  between  100  and  1700 
were  used  to  analyze  the  heat  transfer  characteristics  of  ten 
discrete  flush  mounted  heat  strips  in  an  opposed  flow 
condition.  For  the  initial  part  of  this  section  only  the  odd 
heaters  for  q"=1050  W/m2  will  be  discussed.  However,  it 
should  be  noted  that  the  even  heater  strips  and  other  heat 
flux  settings  exhibited  similar  trends.  Plots  to  be  presented 
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Figure  6.   Natural  convection  comparison  between  Ref[15], 
Ref[16]  and  1181,2  and  3. 


2r> 


are:  Delta  Temp  (THTR-TIN)  vs  Time;  Nu  vs  Re;  and 
Nondimensional  Temp  (©)  vs  Re. 

Figures  7a-l  show  Delta  Temp  vs  Time  for  Reynolds  numbers 
between  1700  and  100  for  odd  heaters  for  q"=1050  W/m2.  These 
figures  show  the  response  of  the  heater  temperatures  with 
decreasing  Reynolds  numbers  for  all  other  parameters  fixed. 
Figures  7a-g  (Re  range  of  1700-500) ,  indicate  conditions  where 
forced  convection  dominates.  At  high  flow  rates  (Figures 
7a, b),  the  temperature  differences  are  collapsed  together  and 
relatively  small  differences  between  the  heat  transfer 
coefficients  at  the  various  locations  exist.  As  flow  rates 
decrease  (Figures  7c-g) ,  a  distinct  separation  between  heater 
strips  is  seen.  The  temperature  difference  is  lower  for  HS9 
(heater  strip  closest  to  inlet  flow)  but  increases  for  HS1 
(heater  strip  furthest  from  inlet  flow) .  Additionally  as  Re 
decreases,  all  AT's  increase. 

For  Re  between  500  and  100,  Figures  7h-j  show  opposed 
mixed  convection  effects — where  buoyant  and  forced  convection 
velocities  are  comparable.  This  is  seen  in  the  almost 
periodic  oscillatory  behavior  of  AT  for  heater  strips  furthest 
from  the  inlet  flow  (Figure  7h)  .  As  Re  decreases,  the 
oscillations  becomes  aperiodic  (Figures  7i  and  j)  with  larger 
amplitudes.  Other  observations  are  the  transposing  of  HS9  and 
HS1  (Figures  7i,j,k)  from  a  minimum  AT  to  maximum  AT  in  the 
heater  array  and  the  more  noticeable  aperiodic  oscillations  of 
HS7  and  HS9  (Figures  7i,j) — the  two  heater  strips  closest  to 

26 


u 

Q. 

E 

<  * 
Q 


01 

6 

2 

0 

4 

0 

fi 

2 

8 

A 

0 


7.CH 

6.6 

6.2-1 

5.8 

5.4 

5.0-1 

4.6 

4  2 

3.8H 

3.4 

3.0 


7.0 
6.6 
6.2 
5.8 
5.4 
50 
4.6 
4.2 
3.8 
3.4 


(a) 


Re  =  1  700 


(b) 


Re  =  1500 


i?TTOro«g^^ 


o 


10 


20 


30 


40 


50 


60 


70 


Mo  -  IL'OO 


30         40         50 

Time  (s) 


HS3 

HS5 
HS7 

X 

HS9 


Figure   7a-c.       Delta    Temp   vs   Time    for    opposed   mixed 
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the  inlet  flow.  It  is  very  clear  that  buoyant  and  forced 
convection  forces  are  competing  for  dominance  between  Re=2  00- 
500.  Figures  7k  and  1  show  transport  dominated  by  buoyancy 
effects  and  the  complete  transposition  of  HS9  and  HS1. 

The  Nusselt  number  variations  for  heat  fluxes  of  1050, 
2500,  3500  and  4500  W/m2  are  shown  in  Figures  8a-d.  As 
observed  in  Figures  7a-g,  a  Re  increase  was  accompanied  by  a 
decrease  in  AT.  For  g"  =  1050  W/m2,  Figure  8a  gives  a  better 
representation  of  this  behavior  in  terms  of  Nu.  For  Re  =  500 
to  1700,  Nu  is  higher  for  heater  strips  closest  to  the  inlet 
flow,  and  it  increases  with  each  heater  strip  as  Re  increases. 
The  ANu  between  HS1  and  HS9  at  Re=500  is  0.95;  at  Re=1700  it 
is  0.17.  This  is  an  85%  reduction  and  shows  a  collapsing  of 
the  data.  Also  note  that  HS1  and  HS9  form  a  high/ low  envelope 
for  all  other  heater  strips.  Another  observation  is  that  as 
Re  decreases,  the  oscillation  becomes  aperiodic  (Figures  7i 
and  j)  and  there  is  a  transition  region  were  Nu  shows  no  clear 
trend  between  Re  =  200  -  500  (Figure  8a) . 

Comparison  of  the  other  heat  fluxes  shows  behavior  similar 
to  Figure  8a.  Heater  strips  9  and  1  envelope  the  other 
heaters  and  there  are  distinct  buoyant,  mixed  and  forced 
convection  dominated  regimes.  Figures  8a-d  show  the 
transition  between  these  three  regimes  and  Table  II  summarizes 
the  mixed  convection  regime. 
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Figure  8a, b.   Nu  vs  Re  for  opposed  mixed  convection  at 
q"=1050  and  2500  w/ra2. 
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Figure    8c, d.       Nu    vs    Re    for    opposed    mixed    convection    at 
q"=3500    and    4500    W/m2. 
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TABLE  II.   MIXED  CONVECTION  REGIME  FOR  OPPOSED  FLOW 


q"  (W/m2) 

Re  range  for  mixed  convection 

1050 

200  -  500 

2500 

300  -  800 

3500 

300  -  1200 

4500 

300  -  1500 

It  is  also  interesting  to  note  the  transposition  of  HS1  and 
HS9  in  Figures  8a-d.  HS1  has  the  highest  Nu  at  low  Re,  but 
the  lowest  at  high  Re.  The  opposite  is  true  for  HS9.  This  is 
reasonable  considering  the  downstream  effects  of  a  row  of 
heaters.  In  buoyancy  dominated  regions,  Nu  for  HS1  will  be 
greater  than  HS9  due  to  the  increasing  thermal  energy 
convected  by  the  flow  resulting  in  higher  temperatures 
downstream  (flow  from  HS1  to  HS9) .  In  a  forced  convection 
dominated  regime,  Nu  for  HS9  will  be  greater  than  HS1  due  to 
the  reversal  in  the  direction  of  the  bulk  flow  (flow  from  HS9 
to  HS1) .   Figures  8a-d  clearly  show  these  trends. 

The  nondimensional  temperatures  (0)  ,  Figures  9a-d,  clearly 
illustrate  the  transition  between  the  transport  regimes  by 
displaying  the  temperature  oscillations  at  HS1  and  HS10  for 
the  four  different  heat  fluxes.  Note  that  the  cross  over 
point  shifts  right  as  heat  flux  increases  and  that  the  heater 
strip  temperature  oscillations  occur  within  the  Re  bands  shown 
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Figure  9a, b.   Nondimensional  Temp  vs  Re  for  opposed  mixed 
convection  at  HS1  and  10  for  q"=1050  and  2500  W/m2. 
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Figure  9c, d.   Nondimensional  Temp  vs  Re  for  opposed  mixed 
convection  at  HS1  and  10  for  q"=3500  and  4500  W/m2. 
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for  q"  =  3500  and  4500  W/m2  do  suggest  trends  similar  to  lower 
q"  data,  but  additional  data  runs  with  Re  >  1700  are  needed 
for  completeness. 

With  the  present  Nusselt  number  data  of  Figures  8a-d,  a 
power  law  curve  fit  was  performed.  Table  III  presents  the 
exponent  (n)  in  the  correlation  Nu  =  Constant*Ren  for  the 
various  heaters. 

Looking  at  the  columns  for  each  heat  flux  shows  the 
collapsing  of  the  data  in  the  forced  convection  dominated 
regime. 

TABLE  III.  EXPONENT  'n'  FOR  THE  POWER  LAW  CURVE  FIT  OF  FIGURES 
8a-d  IN  THE  FORCED  CONVECTION  REGIME 


HEATER 
STRIP 

1050  W/m2 
Gr*=7  0.0k 
Re  ^  500 

2500  W/m2 

Gr*=166.7k 

Re  £  800 

3500  W/m2 
Gr*=2  3  3.4k 
Re  £  1200 

4500  W/m2 

Gr*=300k 

Re  ^  1200 

1 

0.274 

0.182 

0.281 

0.483 

3 

0.229 

0.167 

0.227 

0.373 

5 

0.204 

0.140 

0.180 

0.577 

7 

0.  166 

0.110 

0.132 

0.175 

9 

0.122 

0.130 

0.086 

0.104 
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Figures  lOa-c  illustrate  the  effect  of  position  on  the 
Nusselt  number.  These  figures  show  the  cross  over  point 
shifting  to  the  right  from  HS10  to  HS1.  This  trend  is  similar 
to  that  seen  in  Figures  8a-d.  At  lower  Re  the  flow  is  natural 
convection  dominated,  therefore  as  q"  increases  Nu  increases 
and  this  is  seen  in  Figures  lOa-c.  For  high  Re,  Nu  =  f(Re) 
and  we  see  a  collapsing  of  the  data.  This  collapse  is  less 
noticeable  for  HS1  than  it  is  for  HS10  because  of  the  effect 
of  upstream  heater  strips  on  downstream  strips. 

D.   OPPOSED  MIXED  CONVECTION  LIQUID  CRYSTAL  VISUALIZATION 

The  surface  temperature  pattern  was  visualized  (Figure  11) 
using  a  liquid  crystal  sheet  which  displayed  shades  of  green 
and  blue  in  the  temperature  range  of  3  0  to  3  5  °C. 

Plumes  were  visible  for  q"=3500  and  4500  W/m2  and  evident 
on  all  heaters  for  Re  >  600.  These  plumes  appeared  to  start 
near  the  bottom  heaters  (HS1-4)  at  very  regular  spacings  and 
move  up  along  the  surface.  Near  the  top  heaters,  the  plumes 
were  considerably  wider  and  interaction  effects  were  evident. 
Time-wise  meandering  of  the  plumes  in  the  spanwise  direction 
was  observed. 

At  HS9,  the  plume  showed  no  distinct  streaks,  it  appeared 
as  one  continuous  plume.  Above  HS10,  thermals  appeared.  A 
thermal  is  a  discrete  buoyant  mass  of  fluid  ascending  or 
descending  in  an  ambient  medium  due  to  a  difference  in  density 
[Ref   17:   p.   691].     Distinct   aperiodic   movements   were 
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Figure  11. 


Liquid  crystal  flow  visualization  at 
q"=3500  W/m2  and  Re=400. 
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visualized  as  the  thermals  developed.  The  development 
sequence  started  at  an  arbitrary  location  on  HS10,  grew  into 
a  mushroom  like  form,  necked  and  then  dissipated  away.  This 
sequence  continued  to  repeat  itself  at  various  locations  on 
the  heater  strip. 

As  flow  increased,  the  height  of  the  mushroom  decreased 
and  spanwise  spacing  between  the  individual  plumes  near  the 
bottom  heaters  decreased.  At  higher  flow  rates  where  forced 
convection  dominated,  the  formation  of  thermals  ceased. 
Additionally,  visualization  of  thermal  plumes  faded  away  and 
at  Re  =  1700  none  were  visible.  For  lower  flow  rates  (Re  < 
3  00) ,  the  thermals  would  mingle  together  and  at  the  top  of  the 
test  section  form  into  a  continuous  plume  that  gradually 
dissipated  some  distance  above  HS10. 

The  observed  liquid  crystal  visualization  agreed 
quantitatively  with  the  periodic  and  aperiodic  oscillations  on 
the  Delta  Temp  vs  Time  plots  (Figures  8a-d)  .  At  low  Re, 
buoyancy  forces  dominated  and  AT  showed  no  time-wise 
fluctuations.  As  Re  increased,  instability  ensued  until 
forced  convection  flow  was  established — one  heater  at  time 
starting  with  the  heater  strip  closest  to  the  inlet  flow. 

E.   AIDING  MIXED  CONVECTION  RESULTS 

The  same  four  heat  flux  settings  and  eleven  Re  values  used 
for  opposed  mixed  convection  were  used  to  analyze  aiding  mixed 
convection.   Plots  to  be  presented  and  discussed  are:   Nu  vs 
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Re  for  odd  heaters  and  HS1,  5,  and  10.  The  results  for  aiding 
convection  showed  a  general  increase  in  Nu  with  Re.  In 
Figures  12a-d,  several  trends  are  noticed.  First,  there  is  a 
decrease  in  slope  for  all  heater  strips  as  heat  flux 
increases.  The  largest  decrease  occurs  for  HS1  and  3  (those 
closest  to  the  inlet  flow) .  Table  IV  collects  the  power  law 
curve  fit  exponent  (n)  in  the  relationship:   Nu=Constant*Ren. 

TABLE  IV.   EXPONENT  'n'  FOR  THE  POWER  LAW  CURVE  FIT  OF  FIGURES 
lla-d  IN  THE  FORCED  CONVECTION  REGIME 


HEATER 
STRIP 

1050  W/m2 
Gr*=7  0.0k 
Re  £  800 

2500  W/m2 

Gr*=166.7k 

Re  £  800 

3500  W/m2 

Gr*=2  3  3.4k 

Re  £  800 

4500  W/m2 
Gr*=300k 
Re  ^  800 

1 

0.269 

0.131 

0.117 

0.096 

3 

0.260 

0.139 

0.123 

0.  102 

5 

0.265 

0.153 

0.130 

0.115 

7 

0.254 

0.157 

0.130 

0.108 

9 

0.222 

0.154 

0.128 

0.109 

Table  IV  does  not  show  a  collapse  of  the  data  at  higher  Re; 
but  rather  an  asymptotic  trend  in  the  downstream  direction  for 
each  q" .  Secondly,  at  lower  Re,  Nu  increases  for  all  heater 
strips  as  q"  increases.   This  was  expected  since  in  the 
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Figure  12a, b.   Nu  vs  Re  for  aiding  mixed  convection  for 

q"=1050  and  2500  W/m2. 
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Figure  12c, d.   Nu  vs  Re  for  aiding  mixed  convection  at 

q"=3500  and  4500  W/m2. 
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buoyancy  dominated  region  Nu  =  f(Gr)  and  Gr*  =  f(q").  At 
higher  Re,  in  the  forced  convection  dominated  regime,  the 
change  in  Nu  decreases.  This  effect  is  barely  noticeable  for 
heater  strips  far  downstream  (HS7,9)  at  q"  =  3500  and  4500 
W/m2.  This  suggests  that  Nu  =  f(Re)  and  is  not  dependent  on 
heat  flux.  Finally,  the  downstream  increase  in  heater 
temperatures  results  in  a  decrease  in  the  Nusselt  numbers. 

Figures  13a-c  show  a  buoyancy  dominated  regime  at  low  Re 
and  forced  convection  dominated  regime  at  high  Re.  This 
figure  clearly  shows  the  decrease  in  Nu  downstream,  and  gives 
a  better  visualization  of  reduced  heat  transfer  due  to  the 
cumulative  thermal  boundary  layer  buildup  effect  on  downstream 
heater  strips. 

At  first  glance,  Figure  13a  is  confusing.  At  low  Re,  the 
plot  is  what  was  expected  (ie.  higher  q"  yields  increased  Nu 
due  to  larger  buoyancy  forces)  .  At  higher  Re,  the  Nu  for 
q"=1050  W/m2  exceeds  all  other  heat  fluxes.  Wilks  [Ref  18: 
pp.  743-752]  discusses  the  effect  of  buoyancy  forces  on  the 
boundary  layer  along  a  semi-infinite  constant  heat  flux 
vertical  plate.  At  the  leading  edge,  forced  convection 
dominates  and  buoyancy  effects  are  of  secondary  importance. 
At  positions  further  downstream,  this  trend  is  reversed  and 
buoyancy  forces  dominate  over  forced  convection.  Figure  14 
uses  Wilks'  correlation  for  the  local  Nusselt  number  (Pr  =  1) 
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Figure  13a-c.   Nu  vs  Re  for  aiding  mixed  convection  at 

various  heater  positions. 


46 


Nux  =  ^Re^^e  Q*(e) 

to  suggest  an  explanation  for  the  trend  observed  in  Figure 
13a.  Q*(€),  in  the  equation  above,  is  the  result  of  the  full 
numerical  integration  of  the  boundary  layer  equation  in  Table 
3  of  Reference  [18].  Figure  14  also  shows  that  the  Nusselt 
number  for  q"=1050  W/m2  is  smaller  than  for  q"=4500  W/m2  at 
lower  Re.   At  higher  Re,  this  trend  is  reversed. 
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IV.  CONCLUSIONS 

This  study  reached  the  following  conclusions  on  the 
influence  of  flow  rate,  heat  flux  and  position  on  convective 
heat  transfer  from  discrete  heat  sources  in  a  liquid  filled 
vertical  channel. 

For  forced  convection  dominated  regime  under  both  opposed 
and  aiding  conditions: 

•  As  flow  rates  increased,  Nusselt  numbers  increased. 

•  Heat  transfer  rates  of  downstream  heaters  were  reduced  as 
a  result  of  lying  in  the  thermal  wake  of  upstream  heaters. 
This  was  due  to  the  cumulative  buildup  of  the  thermal 
boundary  layers. 

For  opposed  flow  only,  with  all  other  parameters  fixed, 
and  q"  increased  between  1050  and  4500  W/m2,  the  following 
conclusions  were  drawn.  At  Re<500  where  buoyancy  forces 
dominate,  higher  heat  fluxes  resulted  in  higher  Nu.  For  Re 
between  500  and  1000,  Nu  decreased  and  the  upper  limit  of  the 
mixed  convection  regime  expanded  from  Re=500  to  1500;  the 
lower  limit  showed  a  significantly  smaller  variation  and  was 
in  the  200-300  range.  In  this  range  mixed  convection  regime, 
flow  was  unstable  and  Nu  fluctuated.  At  Re>1000,  a  forced 
convection  regime  was  found  where  Nu  was  a  function  of  Re  and 
independent  of  q" . 

For  aiding  flow,  as  q"  increased  between  1050  and  4500 
W/m2,  Nu  decreased  for  each  heater  strip  and  there  were  no 
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distinct  regimes  as  there  were  for  opposed  flow.  Additionally 
the  slopes  of  the  Nu  vs  Re  curves  decreased  as  q"  increased. 
In  general,  this  study  showed  the  transition  between 
buoyancy  dominated  to  mixed  and  forced  convective  heat 
transfer  to  be  very  unstable  for  opposed  flow,  but  stable  for 
aiding  flow.  Additionally,  for  aiding  flow,  buoyancy  forces 
acted  to  enhance  the  rate  of  heat  transfer  at  low  q" ,  while 
opposed  flow  acted  to  decrease  the  rate  of  heat  transfer. 
This  difference,  however,  decreased  at  higher  heat  fluxes. 
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V.  RECOMMENDATIONS 

In  continuation  of  this  study,  it  is  recommended  that 
follow  on  studies  include: 


•  Velocity  field  determination  using  a  Laser  Doppler 
Velocimeter 

•  Development  of  a  numerical  model  to  predict  the  effects  of 
various  channel  aspect  ratios,  heater  strip  sizes  and 
fluids 
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APPENDIX  A 
UNCERTAINTY  ANALYSIS 


An  uncertainty  analysis  was  performed  to  evaluate  the 
accuracy  of  the  data  in  this  study. 

The  uncertainty  of  a  function  F  =  F(X1,  X2,  X3)  ,  where 


X 


2  3  are  the  independent  measurements  is  determined  by: 


bF  = 


(  Af        \ 
lI  dXi 


dX, 


(  dF 


U/2 


6X3 


V  ^3     /J 


For  F  =  C  X^  X2b  X3C  the  uncertainty  can  be  put  in  the 
following  form: 


M 

F 


bxA2 
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'3f 

X, 


bx3\2 


x, 
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1.   Nusselt  Number  Uncertainty, 


N  = 


hh 


and 
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bN, 


u    _ 


Nr] 


bL\2 
L 


bh\2 
h 


1/2 


where, 


bh 

h 


bQ\2  +  /_6^\2  +  /  SAT 


A  / 


AT 


1/2 


For 


Q   =    2.100    W 


A   =    0.002    m" 


AT    =  3.4    °C 

L   =  0.01    m 

k    =  0.613   W/m  K 

h  =  308.82  W/m  K 

Nu  =  5.04  ±  0.2  6 

The  uncertainty  SNu/Nu  =  0.052  or  5.2%. 


5Q  =  0.0001  W 

6A  =  2.50E-7  m2 

SAT  =  0.05  °C 

6L  =  0.0005  m 


Reynolds  Number  Uncertainty, 


Re=^» 


bR 


E     _ 
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M 


[\      UM 


{    Dh) 
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where, 


bD 


H     _ 


Dg 


6A\2   +  /_6P 
A  /    \    P 


1/2 


For 

A  =  0.002  m2 
P  =  0.42  m 
U  =  0.0767  m/s 

v  =  0.857E-6  m2/s 
Re  =  1700  ±  111 
The  uncertainty  6Re/Re  =  0.065  or  6.5%. 


5A   =  2.5e-7  m2 
6P   =  0.0005  m 


5U  =  0.005  m/s 


Flux  Based  Grashof  Number  Uncertainty, 


Gr*  = 


gVQL4 

KV2A 


bGi* 
Gr* 


f  r + («¥ 


5A  \2 


1/2 


For  the  same  values  given  in  subsection  1  except  6L=0. 000127m 
and  p=306E-6  K"1,  Gr*  =  70,009  ±  3,570. 
The  uncertainty  6Gr*/Gr*  =  0.051  or  5.1%. 
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APPENDIX  B 
SAMPLE  CALCULATIONS 

The  following  calculation  is  for  HS9  where 
q"=1050  W/m2,  Re=1700,  T,N=26.85  °C  and  AT=3  . 4  °C. 

1.  Characteristic  Dimensions. 
Perimeter  (P)  =  2(0.2  +  0.01) 

=  0.42  m 
Heater  surface  area  (As)  =  (0.20) (0.01) 

=  0.002  m2 
Channel  cross  sectional  area  (Ac)  =  (0.02) (0.01) 

=  0.002  m2 
Hydraulic  diameter  (DH)  =  4 (0. 002) /0. 42 

=  0.019  m 

2.  Convective  Heat  Flux. 
Power  to  heater  =  2.100  W 
Heat  flux  (Q")  =  2.100/0.002 

=  1050  W/m2 

3.  Water  Properties  [Ref  8:  p.  A22]. 
P  =  306E-6  K"1 

v  =  .857E-6  m2/s 
k  =  0.613  W/m-K 

4.  Reynolds  Number. 

Re  =  tyyv 

=(0.0767) (0.019) /0.857E-6  =  1700 
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5.  Heat  Transfer  Coefficient. 

h  =  q"/AT 
=  1050/3.4 
=  308.8  W/m2-K 

6.  Nusselt  Number. 
Nu  =  hL/K 

=  308.8(0.01) / (0.613) 
=  5.04 

7.  Fluxed  Based  Grashof   Number. 
Gr*   =   gpq"LA/KV2 

=       (9.81)  (306E-6)  (1050)  (0 . 01) 4/  ( 0  .  613 )  (0.857E-6)2 
=    70,009 
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